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' The Sole-Manrubia model of "macro" -evolution describes the origination, evolution and extinction 

\ of species on geological time scales. We report some properties of this model which provide deep 

. insight into this simple model which captures several realistic features. 
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^ ! I. INTRODUCTION 

' Several fundamental questions on the origination, evolution and extinction of species have been addressed in the 
' recent years by theoretical treatments of dynamic network models using the concepts and techniques of non-equilibrium 
statistical mechanics |^. Kauffman ||^ pioneered a model which describes biological evolution as a "walk" in a 
"fitness landscape" and argued that the eco-system adaptively evolves so as to be delicately poised "at the edge 
^ . of chaos" . On the other hand, Bak and Sneppen 0| , suggested an alternative way of implementing the Darwinian 
' principle of selection (the "survival of the fittest"); this dynamics naturally leads the eco-system to a "self-organized" 
I , critical statejsj. These models implicitly assume a primarily biotic cause of extinctions. However, models based on 
' purely abiotic causes, which arise from "environmental stresses" also lead to similar pattern of extinctions of species 

. ■ A few years ago Sole and Manrubia (SM) ||7[] introduced one such model of "macro" -evolution that takes into 

"j^ ' account the inter-species interactions explicitly. They reported several statistical properties of the model as well as 

: their biological implications. In this paper we report some other statistical properties of this model; some of these 

^ ' features are qualitatively similar to the corresponding results obtained in a more recent model Moreover, we 

'"^ , re-interpret some of the quantities in their approach to extract some relevant biological implications of the model. 

^ Furthermore, we introduce a new correlation function throwing light on the "volatility" of the fitness and the level of 

O , "coherence" in the changes of fitness of the species in this model. 
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II. THE MODEL 



For the sake of completeness, we briefly present the SM model and its main published results. 
The system consists of N species, each labelled by an index i {i = 1,2,...N). The states of the z-th species is 
— [ , represented by a two-state variable S^; Si — or I depending on whether it is extinct or alive, respectively. The 
inter-species interactions are captured by the interaction matrix J; the element denotes the influence of the species 
j on the species i. If Jij > while, simultaneously, Jji < then i is the predator and j is the prey. On the other 
hand, if both Jij and Jji are positive (negative) they the two species cooperate (compete). 

The dynamics of the system consists in updating the states of the system (i.e., to determine the state at the time 
step t + 1 from a complete knowledge of the state at the time t) in the following three steps: 



N 



s,it + l) = e(J2J^sSJit)-o^ (1) 



Step (i): One of the input connections Jij for each species i is picked up randomly and assigned a new value drawn 
(— I from the uniform distribution in the interval [—1, 1], irrespective of its previous magnitude and sign. 
Q Step (ii): The new state of each of the species is decided by the equation 
O 

> 

X. 

^ , where 6i is a threshold parameter for the species i and 8(a;) is the standard step function, i.e., 8(x) = 1 if a; > but 
■ - - • zero otherwise. If S(t + 1) becomes zero for m species, then an extinction of size m is said to have taken place. 

Step (in): All the niches left vacant by the extinct species are refilled by copies of one of the randomly selected 
non-extinct species. 

Sole and Manrubia recorded extinctions of sizes as large as 500 and the distributions of the sizes of these extinctions 
could be fitted to a power law of the form N{m) ^ to^" with an exponent a ~ 2.3. Moreover the periods of stasis 
tg were also found to obey a power law N(ts) ~ tj'' with the exponent 7 ~ 3.0. However, surprisingly, in none of 
their papers H], Sole and his collaborators reported the distributions of the lifetimes of species which, according to 
some claims ^e, for example, refs.|^, |^ for references to the experimental literature and data analysis), also follows 
a power law. 
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FIG. 1: Log-log plots of the distributions of the lifetimes of the species in the SM model with A'^ = 100, ti, 



10^ tr, 



8 X 10® 



(+), N = 300, = tm = 8x 10* (x), N = 1000, t„ 



10^ t„ 



8 X 10 (*). The solid line (with slope —1) and the 



dashed line (with slope —1.5) are drawn just for the comparison of the simulation data with power law of the form Nt ~ t 
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FIG. 2: Distribution of the fitnesses of the species in the SM model with N = 100, tw = 10^ U 
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III. RESULTS AND DISCUSSIONS 



In most of the food webs recorded so far only the sign (but not the magnitude) of the inter-species interactions are 
noted. Moreover, all the qualitative features of the SM model remain unchanged if, instead of allowing the elements 
Jij to take any real value in the domain [—1, 1] with equal probability, we allow to be only +1 or —1 with equal 
probability. Therefore, we work with the simpler version of the SM model where each of the elements of J are discrete 
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FIG. 3: Distributions of tlie incremental changes (per unit time step) in the fitness of the species in the SM model with 



iV = 100, 



10^ t„ 



binary variables that can be either +1 or —1. 

In our simulations we let the system "warm up" for time steps, starting from a random initial condition, and 
then run it for further tm time steps during which collect the data. 



A. Distribution of lifetimes 

It has been claimed in the literature j|, ^ that the distributions of the lifetimes of the extinct species, estimated 
from the fossil data, can be fitted to a power law of the form 

Pt{r)-T-^ (2) 

where several possible values of x, ranging from 1 to 2 have been suggested. However, the log-log plot of our data in 
figj^, obtained from the simulation of the SM model, demonstrates that an approximate power-law-like behaviour is 
seen only over nearly one decade of short lifetimes while strong deviations from power-law are seen for longer lifetimes. 
While the sharp drop of the curves in the tail region can be attributed to finite time effects, the deviations from the 
early power-law-like region at around r ~ 10 is neither an artefact of finite time of simulation nor that of finite size 
of the system |^. 



B. Fitness 



It was originally pointed out by SM Q that fi = J2j Jij^j ~ may be taken as a measure of the "fitness" of the 
species i; the larger the fi the higher is the fitness of the corresponding species i. This definition of "fitness" is, of 
course, consistent with the expectation in the Darwinian scenario that, at any arbitrary stage of evolution, the species 
with positive fitness would not become extinct. Instead of studying the statistics of P{fi), or its time-average, SM 
investigated the relation between the species-averaged value of fi, i.e., fr = J2i fi/-^^ ''^^'^ the size of the extinctions. 

We have now computed the time-averaged probability distribution < P{fi) >t (following SM, we also take 9i — Q 
for all i). From this distribution (see fig.^ it is clear that, on the average, there are very few species which do 
not possess adequate fitness required for survival and, therefore, become extinct. However, the distribution P{fi) 
fluctuates with time and exhibits significant negative fitnesses when extinction occurs. In these stochastic network 
models of ecosystems, the species become extinct because of "bad luck" rather than "bad genes" |10|]. 
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FIG. 4: The sizes of the extinctions and the volatility are plotted against time in (a) and (b), respectively. Data are shown 
only over approximately two and half orders of magnitude, from the top, to highlight the strong correlation between the two 
quantities plotted in (a) and (b). 



We define the incremental change of fitness of a species, say, i, per unit time step as 

AMt) = Mt) - Mt - 1) 

The distribution of Afi is shown in fig.@; it falls roughly exponentially in the tail region. 
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C. Volatility and Coherence 

Borrowing the concept of "volatility" in financial markets |l^, wc now define the volatility" in the eco-system 



V{t)^C,J2{Af,y 



(4) 
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FIG. 5: The correlation function C isplotted against time. The data for C < 10^ are not shown to make the plot less crowded 
and making the comparison with tigm more convenient. 

and the time-dependent correlation function by 

C{t) = C^Y^iAmAf,) (5) 

where the normalization factors Ci and C2 can be chosen appropriately. For simplicity, we worked with Ci = 1 = C2. 
We have also measured the sizes of the extinctions by the fractions m/N where m is the numbers of species that 
become extinct simultaneously at the same time step and N is the total number of species in the system. 

Large extinctions occur when volatility is also large (see fig.§). C is a measure of coherence" in the changes in the 
fitness of the species during the same period of time. Interestingly, this correlation function is very rarely negative. 
Even when it takes negative value, the corresponding magnitude is not large. So, most of the species do not profit 
from the extinction of others. 



D. Single common ancestor 

Each node i {i = 1,2,...,N) carries a label L{i) that indicates the identity of its ancestor. In the initial state we 
assign L(i) — i. During the evolution of the system, whenever a node /3 becomes extinct and is, therefore, re-populated 
by the species occupying the node a we re-label the node /? with L{a), i.e., the node /? also carries the label L{a). 
From our simulation we have always found that, after sufficiently long time, all the nodes carry the same label thereby 
indicating that all the living organisms at that stage share a single common ancestor from which all the wide variety 
of species have appeared as a result of massive diversification. This is perfectly consistent with theory ]l3t and the 
fossil record The distributions of the sizes of phylogenic trees had already been studied by Sole and Manrubia 

0- 

IV. SUMMARY 

We have computed the distribution of the lifetimes of the species in SM model; the deviations from the pure power- 
law behavior observed for the SM model is similar to that in a recent model that describes both "macro" -evolution as 
well as "micro" -evolution ||]. We have also interpreted fi = J2j Jij^j the fitness of the species i and investigated 
the distributions of fitness as well as the incremental changes of fitness per unit time. Finally, we have introduced the 
concepts of volatility and coherence; computing these for the SM model as functions of time, we have demonstrated 
their close relation with the sizes of extinctions. 
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